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We present a swarm model of Brownian particles with harmonic interactions, where the individuals 
undergo canonical active Brownian motion, i.e. each Brownian particle can convert internal energy 
to mechanical energy of motion. We assume the existence of a single global internal energy of the 
system. Numerical simulations show amorphous swarming behavior as well as static configurations. 
Analytic understanding of the system is provided by studying stability properties of equilibria. 



I. CANONICAL ACTIVE BROWNIAN MOTION 



/«~s The concept of active Brownian motion was developed to describe complex motion in various systems [I] , [2] . Its 

^sg main idea is to assume the existence of an additional internal degree of freedom, called the 'internal energy' e, which 

can be taken up from the environment and converted to energy of motion. Former models only assumed an exchange 

' ^ ' between the internal energy and the kinetic energy of motion. In [3], we discussed a more general case, where e can 

Q be also converted to the potential energy of the particle, thus to the full mechanical energy H = ^^ + U{x). The 

^5 corresponding stochastic equations for the position x G R and the momentum p of the particle read 
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1-^ where rj represents the stochastic noise, satisfying the correlations 

§ {vn=0, {v''{t)v\i))^2S'^'Sit-i) (3) 

O 

'— ' and therefore describing a completely irregular force. The constants di and c remain to be fixed for concrete 

, applications. 

> 

\^^ The above system of equations defines the dynamics of a particle subject to an external potential U{x) and 

^^ an irregular force rj. The friction with the environment is characterized by the friction constant 7. The dynamical 

^C equation for e has the following meaning: The first term gives the constant rate of internal energy take-up from the 

CO environment, the second term — ce models the inner energy loss of the indivual, and the third term — ei/describes 

l/~j the conversion of internal energy e into the mechanical energy H. Numerical simulations and analytic investigations 

f^ of these dynamics were made in [T for a single particle in harmonic external potentials. Studies of the multi-particle 

(^^ case were missing so far. 

^ In |3j many particle systems with active Brownian motion were discussed and a specific form of circular swarming 
behavior was found. In this case, however, internal energy was converted only to kinetic energy (corresponding to 

J, di = and H being replaced by ^ in the equation for e). Further simplifications were made by studying only 

C^ situations, where the internal energy e was considered to be stationary (^ = 0) and thus could be eliminated from 
the dynamical equations. 

We now want to discuss in full generality the behavior of many particles undergoing canonical active Brown- 
ian motion, i.e. study particles which can convert dynamical internal energy to total mechanical energy. We will 
show amorphous swarming behavior as well as static configurations. 
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II. MANY PARTICLE SYSTEMS 

We present a simple model for a swarm of active Brownian particles by assuming the existence of a single global 
internal energy e of the swarm. Consider the movement of n active Brownian particles, enumerated by the index 
i, with unit mass in d space dimensions coupled to the internal energy e. The position vectors are x^ G R with 



1, . . . , n and the stochastic equations of motion read 
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The potential energy of one particle is composed of an potential U^ , describing the interaction with all the other 
particles, and an external potential C/j , modeling the environment of the swarm. We study the case where all 
individuals are attracted to the center of mass X = - X]r=i -^i ^^"-^ ^^'^ moving in a landscape of harmonic shape. The 
total mechanical energy therefore reads 

(with positive constants ki and fc2) so that the swarm dynamics is finally given by the equations 

Xj = P^, (9) 

p, = -(l-rfie)(fci(x, -X) + fc2X,) -(7-^26) P, + I7^, (10) 

e=l---eEfY+E^(-^"-.)' + Y-.^V (11) 

2=1 \ j=l J 

A. Variable transformation 

The above system of coupled stochastic nonlinear differential equations is not easily accessible with direct analytic 
procedures. Nonetheless, predictions for the long time behavior can be made by transforming to the new variables 

n 2 ^2 ^ 

^ = Ey' ^-Ey' ^-E-^p^' (12) 
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i^l^Ep^j , ^=(E-«j ' ^=(^E-^j(^Ep^j- (13) 

K represents the total kinetic energy of the swarm, U the total external potential energy, whereas the variables V 
and L are quadratic forms of the center of mass X, resp. its time derivative P = ^: V = ri^X."^ and L ~ n^P^. 



Additionally, also the variable T can be written in terms of X and P: T — n^XP. In these new variables, when 
omitting the noise influences, the system reads: 

k ^ -{1- die) ({ki + k2)S- ^t\ - 2 (7 - d2e) K, (14) 

tl ^ S, (15) 

S = 2K-{1- die) (2 (fci +k2)U- '^v\ - (7 - d2e) S, (16) 

e ^ l-ce-e(K + {ki+k2)U-ki^Y (17) 

T = L-{l^die)k2V-{-f-d2e)T, (18) 

V ^ 2T, (19) 

L = -2{l-die)k2T-2{j-d2e)L. (20) 

B. Equilibrium solutions 

To handle the above nonlinear system, we first calculate stationary points, i.e. values of z = {K,U, S,e,T,V, L), 
where z — and linearize the equations around these equilibrium solutions. By calculating eigenvalues of the Jacobian, 
we can analyze stability properties. We find three distinct equilibrium points Ef. 

El : eo = - , Ko = Uq = Sq = eo = Tq = Vq = Lq = 0. 
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Equilibria Ei and E2 are truly static configurations of the swarm, since Kq = in both cases, whereas equilibirum 
E^ has a nonvanishing equilibrium value of Kq^ which means that swarm particles are actually moving. 

Calculating the stability conditions of all three equilibria allows us to decide under which parameter condi- 
tions the swarm either collapses to a point (Ei), freezes to a certain pattern {E2), or moves amorphously {E3). For 
equilibrium Ei , all seven eigenvalues of the Jacobian are accessible and it is found to be stable under the conditions 

7>— , Odi, 4c(c-di)(fci + fc2) < (7c-c^2)^ (21) 

c 

If these inequalities are satisfied, all swarm particles will collapse into the origin and will remain motionless apart 
from stochastic fluctuations. The evolution of such a collapse, according to equations (14|-(20l is shown in figure IT] 



More interesting is equilibrium E2, which is stable under the conditions 7 > -^, c < di and 

2 {di - c) (fci + fca) < {dij - d2) (^1 + T - J) • (22) 

We used the Routh-Hurwitz-Theorem [5] to determine the stability properties. Exactly one eigenvalue is identical to 
zero independet of the parameter choice. E2 is in fact not an isolated fixpoint, but a one-dimensional continuum of 
equilibria, since Uq and Vq are linearly dependent. This corresponds to exactly one zero-eigenvalue of the Jacobian, 
whereas the other six eigenvalues are stable under the conditions given above. According to stable manifold theory, 
for every chosen value of Vq, there exists a stable manifold on which the trajectory converges asymptotically to 
equilibrium £^2- (In the following, we have set Vq = 0.) In this second equilibrium state, the whole swarm is 
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FIG. 1: Time evolution of the deterministic equations (14l-(20l for K (dashed line), U (solid-line), 5* (dot-dashed line), e 
(solid line), T, V and L (not shown for simplicity of presentation) for the collapse of the swarm to the minimum of the external 
potential (-Bi). Simulation parameters are c — 1, di — 1/2, d2 = 1, 7 = 2, fci = 1/4 and fc2 = 1/4 (Vb = 0). 




- Un 



- "0 



Sr-Ko^O 



FIG. 2: Relaxation to a jelly-like state {E2), according to equations (14l-(20l 
da = 1, 7 = 2, fci = 1/4 and fca = 1/4 (Vb = 0). 



Simulation parameters are c = 1, di = 2, 



frozen to a static configuration, which satisfies Yl^=i ^f = const. One example of such a frozen state is shown 
in figure|2] If the swarm is subject to noise influences, the individual particles will oscillate around their fixed positions. 

Dynamic states of the swarm are reahzed when reaching equihbrium E3, since the sum of all kinetic energies 
then has a non-vanishing value. E^ is found to be stable under the conditions di < 0, 7 < — and 



7di (fci + fcz) < 



4{d2 



di7) {2dl - 2did2j + cdi-f) 



(4d3 - 67^1^2 + 372^2^2 - 073^2) 



(23) 



Again, one zero eigenvalue appears and the Ruth-Hurwitz-Theorem was applied to determine the stability of the 
system. Note that in E3, the two equations for Uq and Lq define two planes in the space of variables {U,V,L), 
whose intersection gives a line, i.e. a manifold of equilibria with dimension one, which matches the overall number 
of zero eigenvalues. Hence the theory of stable manifolds again ensures that we have converging trajectories to E^ 
under the conditions given above. If these inequalities are satisfied by a specific parameter choice, the swarm reaches 
equilibrium E3 and appears to simulate the amorphous behavior of insect swarms. In figure [3J the equilibrating 
process is explicitly shown. 



Numerical simulation of our original equations (l9])-(ll) gives the full iV-particle time evolution. Video exam- 
ples for the corresponding motion of swarm particles are available online [6] for all three distinct cases Ei (collapse) , 
E2 (jelly-like-state) and E^ (amorphous swarming). Simulations of the stochastic system ([9])-(ll) in two dimensions 
are shown, which - apart of small noise perturbations - agree nicely with our simulations of the deterministic 




FIG. 3: State of amorphous swarming (£3), according to equations (14l-(20l. Simulation parameters are c = 1/3, di 
da = 1, 7 = 1, fci = 1/16 and fca = 1/128 (Vb = 0). 



equations (14)-(20) 



III. SUMMARY AND OUTLOOK 

We have formulated swarm dynamics based on canonical active Brownian motion. We assumed the existence of 
a single global internal energy of the swarm and postulated its coupling to the full Hamiltonian of the system. We 
studied the case where all particles were harmonically attracted to their center of mass and were all moving in a 
landscape of harmonic shape. The system was simulated numerically and analytic understanding was obtained by 
studying stability properties of equilibria. We found three distinct forms of the large time behavior of the swarm, 
showing collapse to a point, freezing to a static configuration and amorphous swarming behavior, respectively. 

An interesting application of our present investigation appears to be the study of particle swarms in given 
(genuinely complicated) potentials and the development of a new algorithm for particle swarm optimization theory. 
Furthermore we mention the modeling of crystals and applications in solid state physics. 
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